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1.  Introduction 

Spectacular  improvements  in  jet  noise  prediction  have  been  performed  in  the  last  decade,  due  to 
the  progress  in  numerical  simulations.  Noting  that  hybrid  approaches  based  on  acoustic  analogies 
present  difficulties  related  to  the  modelling  of  source  terms  and  to  the  formulation  of  a wave 
operator  accounting  for  acoustic-flow  interactions,  methods  for  computing  the  sound  directly  from 
the  resolution  of  the  unsteady  compressible  Navier-Stokes  equations  have  been  developed.  The 
objective  of  this  approach  is  to  determine  both  the  flow  field  and  the  acoustic  waves  in  the  same 
calculation.  In  this  way,  the  computed  sound  field  is  a priori  exact  because  no  acoustic  model  is 
used.  It  will  also  permit  to  investigate  the  sound  generation  mechanisms  since  all  flow  and  acoustic 
quantities  are  provided  by  the  computation.  However,  serious  numerical  issues1  must  be  taken 
into  account  in  this  direct  acoustic  approach,  owing  to  the  great  disparity  of  levels  and  length 
scales  between  the  flow  and  the  acoustic  field.  To  overcome  these  difficulties,  numerical  techniques 
specific  to  Computational  AeroAcoustics  (CAA),  suited  to  the  behaviour  of  acoustic  waves,  have 
been  proposed,  such  as  non  dispersive  and  non  dissipative  numerical  schemes,  or  non-reflective 
boundary  conditions.  To  make  direct  aeroacoustic  calculations,  the  challenge  is  then  to  combine 
these  CAA  techniques  with  one  of  the  methods  used  to  solve  the  Navier-Stokes  equations. 

The  three  classical  approaches  commonly  used  to  solve  these  flow  equations  have  thus  been 
experimented  for  aeroacoustic  simulations:  Direct.  Numerical  Simulations  (DNS)  consisting  in  com- 
puting .all  turbulent  scales,  Large  Eddy  Simulations  (LES)  where  only  larger  scales  are  calculated 
whereas  the  effects  of  unresolved  ones  are  modelized  via  a subgrid  scale  model,  and  the  unsteady 
Reynolds  Averaged  Navier-Stokes  equations  (RANS)  using  turbulence  closures.  Among  the  first 
three-dimensional  applications  performed  in  the  last  five  years,  we  can  put  forward  the  DNS  of 
Freund  et  al.2  :i  of  the  noise  radiated  by  supersonic  and  subsonic  circular  jets,  the  LES  of  Morris  et 
a/.4  of  the  radiation  of  a supersonic  rectangular  jet,  and  the  study  of  screech  tones  generation  in  a 
round  jet  using  unsteady  RANS  by  Shen  & Tam.5 


Paper  presented  at  the  RTO  AVT  Symposium  on  "Ageing  Mechanisms  and  Control: 
Part  A -Developments  in  Computational  Aero-  and  Hydro-Acoustics", 
held  in  Manchester,  UK,  8-11  October  2001,  and  published  in  RTO-MP-079(I). 
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The  use  of  LES  is  important  because  it  can  theoretically  be  applied  to  high  Reynolds  number 
flows  unlike  DNS  which  is  restricted  to  very  low  Reynolds  number  flows.  Therefore,  in  the  present 
study,  a Large  Eddy  Simulation  of  a three-dimensional  circular  jet  with  a Mach  number  of  0.9  and 
a Reynolds  number  of  65000  is  performed,  carrying  on  earlier  preliminary  works.6  The  motivation 
is  to  show  the  feasibility  of  the  the  direct  computation  by  LES  of  the  sound  generated  by  a sub- 
sonic flow.  For  this,  CAA  numerical  methods  are  used  in  the  simulation.  The  validation  of  the 
computation  is  then  carried  out  by  comparing  systematically  both  flow  properties  and  the  radiated 
sound  field  with  experimental  data. 

This  paper  is  organized  as  follows.  Numerical  procedure,  characteristics  of  the  jet  and  parame- 
ters of  the  simulation  are  presented  in  section  2.  The  flow  development  is  briefly  shown  in  section 
3,  by  studying  the  mean  flow  and  the  turbulent  intensities.  The  radiated  acoustic  field  is  then 
investigated  hi  section  4.  Finally,  concluding  remarks  are  given  in  section  5. 

2.  Numerical  simulation 

2.1  Numerical  procedure 

The  full  three-dimensional  Navier-Stokes  equations  written  in  a conservative  form  are  solved. 
The  numerical  algorithm  is  low-dispersive  and  low-dissipative  to  compute  directly  the  sound  waves 
with  accuracy,  and  has  been  successfully  used  in  an  earlier  study  on  a two-dimensional  mixing  layer.6 
It  combines  the  Dispersion-Relation-Preserving  (DRP)  finite-difference  scheme  of  Tam  & Webb7 
for  space  discretization  with  a four-step  Runge-Kutta  scheme  for  time  integration.  A selective 
damping1  is  also  applied  to  filter  out  numerical  oscillations,  typically  grid-to-grid  waves,  which  are 
not  computed  properly  by  the  algorithm,  and  to  ensure  numerical  stability.  The  mesh  is  Cartesian 
and  non  uniform,  to  use  different  discretizations  in  the  flow  field  and  in  the  acoustic  far-field. 

The  resolution  of  the  full  Navier-Stokes  equations  without  modelling,  by  Direct  Numerical 
Simulations  (DNS),  requires  the  calculation  of  all  turbulent  scales  and  is  intrinsically  restricted  to 
low  Reynolds  number  flows.  To  simulate  flows  at  higher  Reynolds  number,  Large  Eddy  Simulation 
(LES)  can  be  performed  by  computing  only  the  larger  scales  and  by  taking  into  account  the  effects 
of  unresolved  ones  through  a subgrid  scale  model.  In  the  present  study,  a turbulent  viscosity  is 
classically  implemented  to  provide  the  dissipation  of  the  unresolved  scales.  Various  models  have 
been  developed  to  give  an  evaluation  of  the  turbulent  viscosity  from  the  resolved  scales.  To  keep 
the  problem  as  simple  as  possible  for  aerodynamics,  we  choose  the  Smagorinsky  model.8 

2.2  Boundary  Conditions 

In  acoustic  simulations,  boundary  conditions  must,  minimize  the  magnitude  of  acoustic  waves 
reflected  when  fluctuations  leave  out  the  computational  domain.  In  this  study,  a formulation9 
derived  from  the  behaviour  of  sound  waves  in  the  acoustic  far-field  is  used.  It  is  the  extension 
of  Tam  & Dong10  conditions  for  a three-dimensional  geometry.  Two  conditions  are  applied:  a 
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radiation  condition  for  the  inflow  and  for  the  lateral  boundaries  of  the  computational  domain  where 
only  acoustic  fluctuations  propagates,  and  an  outflow  condition  for  the  outflow  where  vortical  and 
entropic  disturbances  are  also  convected  by  the  flow. 

The  efficiency  of  the  boundary  conditions  has  been  estimated  using  test  cases.13  The  exit  of 
vortical  structures  generates  low-amplitude  spurious  waves,  which  are  however  not  negligible  with 
respect  to  the  physical  sound  field  characterized  by  very  small  amplitudes.  Therefore  a sponge  zone 
is  implemented  in  the  outflow  direction  to  dissipate  flow  fluctuations  before  they  reach  the  outflow 
boundary,  and  also  to  filter  out  possible  reflected  waves.  It  is  based  on  the  combination  of  grid 
stretching  with  an  artificial  damping/’ 

2.3  Flow  Parameters 

The  inflow  longitudinal  velocity  u (r)  of  the  jet  is  given  by  the  following  hyperbolic-tangent 
profile 

where  Uj  is  the  inflow  centerline  velocity,  Sg  the  initial  momentum  thickness  of  the  shear  layer,  and 
7'o  the  jet  radius.  The  fluid  surrounding  the  jet  is  initially  at  rest. 

The  jet  Mach  number  is  Mj  = 0.9.  This  choice  is  justified  by  the  amount  of  experimental 
studies  available  in  the  literature,  providing  both  aerodynamic  results,12  and  acoustic  results. 1:1-15 
This  Mach  number  implies  also  a high  convection  speed  of  turbulent  structures,  which  reduces 
computation  time.  We  can  notice  that  the  first  simulation  of  a subsonic  jet  to  determine  directly 
its  radiated  field,  a DNS  carried  out  by  Freund,3  involves  a Mach  number  0.9,  Reynolds  number 
3600  jet.  In  the  present  simulation,  the  jet  Reynolds  number  is  ile©  = Uj  x D/v  — 65000.  It  is 
higher  than  Reynolds  numbers  reachable  by  DNS,  but  still  lower  than  values  of  practical  interest. 
The  ratio  between  the  initial  jet  radius  and  the  initial  momentum  thickness  of  the  shear  layer  Sq/tq 
is  also  an  important  parameter,  because  the  laminar-turbulent  transition  depends  on  it.  In  this 
study,  dg/ro  = 0.05,  and  the  development  of  vortical  structures  in  the  shear  zones  is  possible. 

Finally,  to  obtain  the  natural  flow  development,  the  jet  is  forced  by  adding  random  velocity 
fluctuations  into  the  inflow.  These  fluctuations  are  solenoidal  to  minimize  the  generation  of  spurious 
acoustic  waves  by  the  excitation.  They  are  introduced  only  in  the  shear  layers,  which  are  the 
unstable  region  of  the  jet.  In  these  shear  zones,  turbulent  intensities  generated  by  the  forcing  are 
around  3%,  which  is  similar  to  intensities  measured  near  jet  nozzle  exits. 

2.'4  Numerical  specifications 

The  computational  domain,  discretized  by  a grid  of  255  x 187  x 127  points,  is  defined  by 
0 < x/ro  < 30,  —5  < y/ro  < 25  and  —5  < z/r o < 5 in  the  three  coordinate  directions.  Radially, 
points  are  clustered  in  the  jet  with  26  points  in  the  radius.  The  mesh  spacing  is  minimum  around 
r — r o with  Aq  = 5g/\.§,  and  increases  outside  the  jet  to  reach  Aymax  = 0.4ro  in  far-field.  In  the 
flow  direction,  the  mesh  spacing  is  constant  up  to  x — 20ro,  with  Ax  = 3Aq.  Then,  meshes  are 
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stretched,  to  obtain  Axmax  = 0.54ro,  to  form  a sponge  zone.  An  artificial  damping  is  moreover 
applied  in  this  zone,  in  the  flow  region  only.  In  this  way,  the  physical  part  of  the  computational 
domain  extends  up  to  x = 20?o  for  the  flow  field,  and  up  to  x = 30?o  for  the  acoustic  field. 

The  time  step  is  given  by  At  = 0.7Ao/co-  The  simulation  runs  for  30000  iterations,  the 
calculation  of  statistical  means  being  performed  only  after  5000  iterations.  Consequently,  the 
simulation  time  corresponds  to  22  x Lx/co , where  Lx  = 30?-o-  The  computation  is  15  hours  long 
on  a Nec  SX-5,  with  a CPU  time  of  0.3/rs  per  grid  point  and  per  iteration.  More  details  of  the 
simulation  can  be  found  in  earlier  papers.9,11 

3.  Flow  Field 


3.1  Flow  Development 


Figure  1 displays  the  vorticity  fields  provided  by  the  simulation.  The  longitudinal  vorticity  field 
uxy  shows  that,  initially,  vortical  structures  are  generated  and  grow  in  the  shear  layers.  The  shear 
zones  of  the  jet  interact  around  x ~ 10ro,  indicating  the  end  of  the  potential  core,  which  is  consistent 
with  experimental  observations.12  A developed  turbulence  is  found  downstream,  illustrated  by  the 
transverse  vorticity  field  ojyz , and  characterized  by  a typical  three-dimensional  mixing.  Finally, 
flow  fluctuations  are  dissipated  by  the  sponge  zone  for  x > 20o- 


y'ro 


Figure  1:  Snapshots  of  the  vorticity  field.  Upper  picture,  u)xy  in  the  x — y plane  at  z = 0;  bottom 
picture,  ioy.  in  the  y — z plane  at  x = 12r0.  Levels  are  given  in  s-1. 
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3.2  Mean  Flow  Properties 

The  characteristics  of  the  jet  mean  flow  calculated  from  the  simulation  are  presented  in  Figures 
2.  First,  Figure  2(a)  displays  the  longitudinal  evolution  of  the  inverse  of  the  mean  centerline 
velocity  normalized  by  the  inflow  velocity  Uj.  Its  value  is  1 in  the  potential  core,  and  grows  linearly 
afterwards  so  that  we  can  write  Uc/Uj  = B x D/(x  — xq)  where  xq  = 0 is  the  virtual  origin,  and 
B = 5.5.  This  value  is  consistent  with  measurements  of  Wygnanski  & Fiedler,10  Panchapakesan  & 
Lumley17  and  Hussein  et  a/.,18  as  well  as  DNS  result  of  Boersma  et  a/.,19  reported  in  Table  1. 


RcD 

B 

A 

Reference 

8.6  x 104 

5.4 

0.086 

Wygnanski*  et  al.u> 

1.1  x 104 

6.1 

0.096 

Panchapakesan*  et  al.17 

"9.5  x 104 

5.8 

0.094 

Hussein*  et  al.Lii 

2.4  x 103 

5.9 

0.095 

Boersma**  et  al.1® 

6.5  x 104 

5.5 

0.096 

Present  simulation 

Table  1:  Mean  flow  parameters  obtained  from  experiments*,  DNS**,  and  present  simulation. 

Second,  Figure  2(b)  presents  the  longitudinal  evolution  of  the  jet  half-width  8l/2  normalized 
with  the  radius  ro.  Its  value  is  fairly  1 up  for  x < 6ro  where  no  large  vortical  structure  is  visible. 
Afterwards,  in  the  turbulent  region,  the  jet  spreads  linearly  as  81/2  = A x (x  — a-'o)  where  A = 0.096, 
that  is  in  agreement  with  data  of  Table  1. 

Third,  Figure  2(c)  shows  the  longitudinal  evolution  of  the  flow  rate.  It  grows  regularly,  indi- 
cating that  the  entrainment  of  the  surrounding  fluid  occurs  from  the  inflow  of  the  computational 
domain.  As  deduced  from  the  a:-1  decay  of  the  centerline  velocity  and  from  the  linear  spreading,  a 
linear  growth  is  found  after  the  potential  core,  as  Q/Qo  = C x (a;  — Xo)/D  where  Qq  is  the  initial 
flow  rate  and  C = 0.32.  This  corresponds  exactly  to  the  measurement  of  Ricou  & Spalding.20 

Finally,  radial  profiles  of  the  mean  longitudinal  velocity  normalized  by  the  centerline  velocity 
Uc  are  plotted  in  Figure  2(d)  as  function  of  the  nondimensional  coordinate  y/(x  — a:o),  for  five 
locations  between  x = 15ro  and  x = 20?'o.  The  spreading  rate  A is  the  half-width  of  the  profiles 
which  are  superimposed.  This  confirms  the  self-similarity  in  the  turbulent  jet  suggested  by  the 
previous  longitudinal  evolutions. 

3.3  Turbulent  Intensities 

Intensities  of  flow  perturbations  in  the  jet  are  investigated.  They  are  calculated  in  the  x — y 
plane  at  z = 0 using  velocity  fluctuations  u',  v'  et  w'  and  normalized  with  the  local  mean  centerline 
longitudinal  velocity  Uc.  They  are  given  by  the  following  expressions:  auu  = vV/(/c,  crvv  — 
y&fUc,  &ww  = \fw^/Uc,  and  auv  = \J\u'v'\/Uc. 

Radial  profiles  of  the  turbulent  intensities  are  shown  in  Figures  3 for  five  locations  between 
x = 15ro  and  x = 20ro,  as  function  of  the  nondimensional  coordinate  y/(x  — #o).  Profiles  are 
well  superimposed,  that  is  in  agreement  with  the  self-similarity  in  a fully  turbulent  jet.  The  mean 
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Figure  2:  Longitudinal  evolution  of:  (a),  the  inverse  of  the  mean  centerline  velocity  normalized 
with  the  inflow  velocity  Uj/Uc\  (b),  the  half-width  of  the  jet  normalized  with  the  jet  radius  <5i2/ro; 
(c),  the  mean  flow  rate  normalized  with  the  inflow  rate  Q/Qo-  (d),  Radial  profiles  of  the  mean 
longitudinal  velocity  normalized  by  the  local  centerline  velocity  U/Uc:  x,  at  x = 15.8r0;  +,  at 

x = l6.8ro;  *,  at  x — 17.7ro;  at  * = 18.6ro;  A,  at  x — 19.67-o.  , approximation  with  a 

Gaussian  profile. 
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profiles  calculated  between  x = 15ro  and  x = 20ro  are  also  plotted  as  solid  lines.  They  are  very 
close  to  measurements,  and  stand  between  two  experimental  profiles  provided  by  Panchapakesan 
& Lumley17  and  by  Hussein  et  al.18  respectively. 

(a)  (b) 


Figure  3:  Radial  profiles  of  the  turbulent  intensities:  (a),  <juu;  (b)  5 &VV  ; (c),  crww.  x,  at  x = 15.8ro; 

+,  at  x = 16.8ro;  *,  at  x = 17.7ro;  V>  at  x = 18.6ro;  A,  at  x = 19.6ro.  , mean  profiles 

calculated  from  x = 15ro  to  x — 20ro; , experimental- profiles  obtained  by  Panchapakesan 

& Lumley17j , experimental  profiles  obtained  by  Hussein,  Capp  & George.18 
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4.  Sound  Radiation 


4.1  Dilatation  field 

Figure  4 displays  the  dilatation  field  © = V.u  provided  directly  by  the  simulation.  Outside  the 
flow  region  where  the  mean  flow  is  negligible,  dilatation  is  proportional  to  the  time  derivative  of 
the  acoustic  pressure  via  the  relation  © = — (l/po^o)  dp/dt. 


Figure  4:  Snapshot  of  the  dilatation  field  © = V.u  in  the  acoustic  region,  and  of  the  vorticity  field 
uxy  in  the  aerodynamic  region,  in  the  x — y plane  at  2 = 0.  The  dilatation  color  scale  is  defined 
for  levels  from  -90  a 90  s-1,  the  vorticity  scale  is  the  same  as  in  Figure  1. 

This  figure  demonstrates  that  the  dilatation  field  is  not  contaminated  by  spurious  waves  gen- 
erated by  the  inflow  excitation  or  reflected  by  the  boundary  conditions.  Moreover,  wave  fronts 
are  mainly  coming  from  an  origin  located  around  x — 10ro,  in  the  region  where  the  mixing  layers 
are  merging,  showing  that  predominant  sound  sources  are  located  in  the  vicinity  of  the  end  of 
the  potential  core.  This  agrees  both  with  the  results  of  the  recent  DNS  performed  by  Freund3 
and  with  the  measurements  of  Chu  & Kaplan21  and  Juve  et  al .22  using  various  source  localization 
techniques.  The  computed  radiated  field  is  also  more  pronounced  in  the  downstream  direction 
according  to  experimental  directivities. 

4.2  Pressure  field  spectra 

To  characterize  the  computed  acoustic  field,  sound  pressure  spectra  are  calculated  in  far  field. 
They  are  transposed  at  a distance  of  60ro  from  the  inflow  using  the  r~l  decay  of  acoustic  waves 
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from  the  sound  sources  located  at  x = 10ro  to  the  recording  points. 

Figures  5(a)  and  (b)  show  the  spectra  found  for  angles  of  30°  and  of  90°  from  the  jet  axis  for 
Strouhal  numbers  St  > 0.15.  Levels  are  more  than  10  dB  higher  for  0 = 30°  than  for  0 = 90°. 
Shapes  are  also  different:  for  0 = 30°,  the  spectrum  is  dominated  by  low  frequency  components 
with  a peak  found  around  St  ~ 0.2,  whereas  for  0 = 90°,  the  spectrum  is  larger  with  a preferred 
band  of  0.2  < St  < 0.5  and  a peak  evaluated  around  St  ~ 0.3.  This  behaviour  corresponds  to 
experimental  results,14,22  and  the  classical  interpretation  is  that,  at  low  angles,  the  noise  is  mainly 
generated  by  the  large  scales,  whereas  at  higher  angles  the  fine-scale  turbulence  is  responsible  for 
a great  part  of  the  noise. 

(a)  (b) 


Figure  5:  Sound  pressure  level  as  function  of  Strouhal  number  St  — fD/Uj,  at  60ro  from  the  jet 
nozzle,  for  an  angle  of:  (a),  0 = 30°:  (b),  6 = 90°. 

4.3  Sound  pressure  level 

The  pressure  spectra  are  now  integrated  from  the  Strouhal  number  St  = 0.15  to  provide  the 
sound  directivity.  Figure  6 shows  the  computed  sound  pressure  levels,  with  experimental  data 
obtained  by  Mollo-Christensen  et  a/.,13  Lush,14  and  Stromberg  et  al.15  for  jets  with  similar  Mach 
numbers  but  various  Reynolds  numbers 

The  agreement  between  calculated  and  measured  sound  levels  is  very  good  for  all  observation 
angles. _ As  expected,  the  acoustic  level  reaches  a peak  around  an  angle  of  0 = 30°.  By  considering 
the  similarities  between  the  noise  radiated  by  jets  with  very  different  Reynolds  numbers  from  3600 
up  to  5.4  x 105,  the  sound  sources  associated  to  the  large  scales  seem  to  be  relatively  independent 
of  the  Reynolds  number.  Moreover,  at  high  observation  angles,  the  result  of  the  simulation  stands 
between  the  measurements  which  are  relatively  scattered.  This  behaviour  can  be  associated  to  a 
Reynolds  number  effect,  because  the  fine-scale  turbulence  is  much  less  developed  at  low  Reynolds 
numbers. 
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Figure  6:  Sound  pressure  level  as  function  of  angle  6 measured  from  the  jet  axis,  at  60ro  from  the 
jet  nozzle.  Experimental  data  by:  +,  Mollo-Christensen  et  a/.13  ( Mj  — 0.9,  Rep  = 5.4  x 105);  o, 
Lush14  (Mj  = 0.88,  Rep  = 5 x 10s);  x,  Stromberg  et  a/.15  (Mj  = 0.9,  Rep  — 3600). 

5.  Conclusion 

Results  of  a Large  Eddy  Simulation  of  a circular  jet  with  a Mach  number  of  0.9  and  a Reynolds 
number  of  65000,  carried  out  to  compute  directly  the  radiated  sound  field,  has  been  presented.  Both 
the  flow  field  and  the  sound  field  computed  have  been  investigated  and  compared  systematically  to 
experimental  data  available  in  the  literature.  Flow  properties,  namely  flow  development,  mean  flow 
parameters  and  turbulent  intensities,  are  in  very  good  agreement  with  measurements.  The  sound 
field  provided  by  LES  compares  also  successfully  with  experimental  results  in  terms  of  directivity, 
spectra  and  levels.  Sound  sources  in  the  jet  are  moreover  found  around  the  end  of  the  potential 
core,  as  shown  experimentally. 

This  study  shows  the  feasibility  of  the  direct  calculation  of  the  acoustic  field  generated  by 
subsonic  flows  using  LES.  The  excellent  concordance  with  measurements  supports  that  aerodynamic 
and  acoustic  mechanisms  are  well  accounted  for  by  the  simulation.  In  that  way,  further  works  would 
have  to  investigate  the  generation  mechanisms  of  jet  noise. 
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